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Intrinsic energy cut-off in diffusive shock acceleration: 
possible reason for non-detection of TeV-protons in SNRs 
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Abstract. The linear theory of shock acceleration predicts 
the maximum particle energy to be limited only by the ac- 
celeration time and the size of the shock. We study the com- 
bined effect of acceleration nonlinearity (shock modification 
by accelerated particles, that must be present in strong astro- 
physical shocks) and propagation of Alfven waves that are 
responsible for particle confinement to the shock front. We 
show that wave refraction to larger wave numbers in the non- 
linearly modified flow causes enhanced losses of particles in 
the momentum range p ma x/ R < P < Pmax, where R > 1 
is the nonlinear pre-compression of the flow and p max is a 
conventional maximum momentum, that could be reached if 
there was no refraction. 



1 Introduction 

One of the most important parameters of the Fermi acceler- 
ation is the rate at which it operates. Indeed, what is often 
predicted or even observed is a power-law spectrum that cuts 
off due to the finite acceleration time. The cut-off momentum 
Pmax (t) advances with time as follows 
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while the acceleration time is determined by (e.g., Axford, 
1981) 
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with m and u 2 being the upstream and downstream flow 
speeds in the shock frame whereas k\ and k 2 are the par- 
ticle diffusivities in the respective media. These are the most 
sensitive quantities here which are determined by the rate at 
which particles are pitch angle scattered by the Alfven tur- 
bulence. If the latter was just background turbulence in the 



interstellar medium, the acceleration process would be too 
slow to account for the galactic cosmic rays (CRs). How- 
ever it was realized (e.g., Bell, 1978) that accelerated par- 
ticles must create the scattering environment by themselves 
generating Alfven waves via the cyclotron instability. This 
wave generation process proved to be very efficient (e.g., 
Volk et al, 1984) so that the normalized wave energy den- 
sity iSB/Bo) is related to the partial pressure P c of CRs 
that resonantly drive these waves through 



(8B/B ) 2 ~ M A P c /pu 2 



(3) 
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where Ma 3> 1 is the Alfven Mach number and pu 2 is the 
shock ram pressure. Often it is assumed that the turbulence 
saturates at SB /Bo ~ 1, which means that the m.f.p. of pitch 
angle scattered particles is of the order of their gyro-radius 
r g . Then, k = kb = cr g (p)/3, where kb stands for the 
Bohm diffusion coefficient. Hence, t acc ~ (eB/p) (c/uj.) . 

However the acceleration rate (||) with k = kb was found 
to be fast enough to explain the acceleration of CRs in SNRs 
up to the "knee" energy ~ 10 15 ey over their life time. The 
analyses of Drury etal, (1994) and Naito and Takahara (1994) 
of prospective detection of super- TeV emission from nearby 
SNR (they must result from the decays of 7r° mesons born in 
collisions of shock accelerated CR protons with the nuclei of 
interstellar gas) look equally optimistic. The expected fluxes 
were shown to be strong enough to be detected by the imag- 
ing Cherenkov telescopes. Moreover, the EGRET (Esposito 
et al, 1996) detected a lower energy (;$ TeV) emission co- 
inciding with some galactic SNRs. Unfortunately, despite 
the physical robustness of the above-mentioned predictions 
of emission, no statistically significant signal that could be 
attributed to any of the EGRET sources was detected (Buck- 
ley et al, 1997). The further complication is that the region 
between GeV and TeV energy bands is currently uncovered 
by any instrument. Therefore, based on these observational 
results it was suggested (e.g., Buckley et al, 1997) that there 
is probably a spectral break or even cutoff somewhere within 
this band. 



In this paper we attempt to understand what may happen 
to the spectrum provided that the acceleration is indeed fast 
enough to access the TeV energies over the life time of SNRs 
in question. Our starting point is that the fast acceleration 
also means that the pressure of accelerated particles becomes 
significant relatively early and must change the entire shock 
structure by this time. At the first glance this should not slow 
down acceleration since according to eq.(g) this increases the 
turbulence level improving thus particle confinement near the 
shock front and making acceleration faster. Simultaneously 
with that but more importantly, the upstream flow is decel- 
erated by the pressure of CRs P c which influences the spec- 
tral properties of the turbulence by affecting the propagation 
and excitation of the Alfven waves. This effect is twofold. 
First the waves are compressed in the converging plasma flow 
upstream and are thus blue-shifted lacking the long waves 
needed to keep the high energy particles diffusively bound to 
the accelerator. Second, and as a result of the first, at highest 
energies there remain relatively few particles so that the level 
of resonant waves is also small and hence the acceleration 
rate is low. 



2 Basic Equations and Approximations 

We use the standard diffusion-convection equation for de- 
scribing the transport of high energy particles (CRs) near a 
CR modified shock 
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then gradually increases up to u\ = u(+oo) > uq (see 
Fig. la). Limiting our consideration to high Mach number 
shocks, M ^> 1, we may drop P g term upstream, x > 0. It 
is retained at the subshock which, however, can be described 
by the conventional Rankine-Hugoniot jump condition 
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where Mq is the Mach number in front of the subshock. For 
simplicity we use the adiabatic approximation, i.e., the far 
upstream Mach number is related to Mq by Mq — M 2 /R y+1 , 
where R = u\/uq is the flow precompression in the CR pre- 
cursor. 

For determining the CR diffusion coefficient n one needs 
to write the wave kinetic equation for which we use the eikonal 
approximation 
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Here jV^ is the number of wave quanta, u> is the wave fre- 
quency u> = —ku + kVA — —ku, k is the wave number. 
The left hand side has a usual Hamiltonian form that states 
the conservation of N k along the lines of constant frequency 
u){k,x) = const on the k,x plane. The first term on the 
r.h.s. describes the wave generation on the cyclotron insta- 
bility of a slightly anisotropic particle distribution. It can be 
expressed through its spatial gradient. The resonance condi- 
tion kp = const ("resonance sharpening," e.g., Drury et at, 
1996 [D96]) is implied. 



Here x is directed along the shock normal (also the direction 
of the ambient magnetic field). The two quantities that con- 
trol the acceleration process are the flow profile U(x) and 
the particle diffusivity K,(x,p). The first one is coupled to 
the particle distribution / through the equations of mass and 
momentum conservation 
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3 Outline of the Analysis 

It is convenient to use the wave energy density normalized to 
dink and to the energy density of the background magnetic 
field -Bq/87t instead of N k 
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along with the partial pressure of CRs normalized to <ilnp 
and to the shock ram pressure p\u\ 
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Using this variables, denoting g = P/p, assuming a steady 
state andp ^> 1, eqs.(Ep) rewrite 



is the pressure of the CR gas, and P g is the gas pressure. The 
lower boundary po in the momentum space separates CRs 
from the thermal plasma and enters the equations through the 
magnitude of / at p = po which specifies the injection rate. 
The particle momentum p is normalized to mc. We assume 
that the upstream region is x > half-space, and represent 
the velocity profile in the shock frame as U(x) — —u(x) 
where the (positive) flow speed u(x) jumps from u^ = u(0— ) 
downstream to uq = u(0+) > 112 across the subshock and 



dx 
dl_ 
dx 



ug 



u x p- 



dx 
d I 



1 dg 
- 3 U * P d~p 
2u\ d_ p 



dpp 2 Va dx 



(12) 
(13) 



Here the CR diffusion coefficient k can be expressed through 
the wave intensity by k = kb/I ■ The difference between 
these equations and those used in e.g., D96 is due to the 



terms with u x ^ 0. Far away from the subshock where u x — > 
0, one simply obtains I = 2u\P/Va. The most important 
change to the acceleration process comes from the terms with 
u x 7^ 0. Indeed, let us recall first how the equation (O) may 
be treated in the linear case u x = for x > 0. Integrating 
both sides between some x > and x — oo, one obtains 
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where we denoted ko = i^b/p — const forp ^> 1. Although 
this equation has a formal spatial scale I ~ Ko/uxM\g, its 
only solution is a power law in x 
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and thus has no scale. It simply states the balance between 
the diffusive flux of particles escaping upstream (second term 
in eq.Jl4|]) and their advection with thermal plasma in the 
downstream direction (the first term). As we shall see, this 
balance is possible not everywhere upstream and the physi- 
cal reason why it appears to be so robust in the case u x = 
is that the flows of particles and waves on the x,p-plane 
(including the diffusive particle transport) are both directed 
along the x-axis. If, however, the flow modification upstream 
is significant (u x > 0, x > 0), the situation changes funda- 
mentally. Fig. 1 explains how the flows of particles and waves 
on the x,p-plane become misaligned even though they are 
both advected with the thermal plasma. In fact, the flows 
separate from each other and, since neither of them can exist 
without the other (waves are generated by particles that, in 
turn, are trapped in the shock precursor by the waves) they 
both disappear in that part of the phase space where the sep- 
aration is strong enough. To understand how this happens it 
is useful to rewrite eqs.(0|Tl) in the following characteristic 
form (we return to the particle number density /) 



u(x) 



(a) 



p 

max 


k P (b) * _^ — 
\ ^^ ^^^^^^ 






/ \<^^^^^^^ ~~~-~- ??£ticles 






/ s' ^ Front 4(x,p)=0 




p, 








"0 


X 









Fig. 1. (a) flow velocity; (b) characteristics of eqs.(|16|jl7|) 



of this paper). Then, integrating eq.(|l2J) between some x > 
and x = oo, instead of (hMb we obtain 
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In contrast to the solution of eq.fjlj) the length scale L = 
Ko/2uiMa enters the solution of this equation. This is be- 
cause it has a nonzero r.h.s. The solution of this equation 
changes rapidly on a scale L -C L c where L c = K.(p max ) / U\ 
is the total scale height of the CR precursor. In addition to x 
and p, we introduce a fast (internal) variable £(x,p) as fol- 
lows 
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One sees from the l.h.s.'s of these equations that particles 
are transported towards the subshock in x and upwards in 
p along the family of characteristics up 3 = const, whereas 
waves move also towards the subshock but downwards in p 
along the characteristics ujp = const. As long as u(x) does 
not significantly changes the waves and particles propagate 
together (along a;-axis) as e.g., in the case of unmodified 
shock or far away from the subshock where u x — > 0. When 
the flow compression becomes important (u x ^ 0) their sep- 
aration leads to decrease of both the particle and wave energy 
densities towards the subshock. To describe this mathemati- 
cally, let us assume that the linear relation between P and / 
is still a reasonable approximation even if u x is nonzero but 
small. (A more general case is considered in a longer version 
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where x = X{(p) is some special curve on the x,p plane 
which bounds the solution and will be specified later. We 
rewrite eq.(|18|) for £ — fixed, L — > which leads to the 
solution 
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where w(p) = Ut + (l/3)pdui/dp, U{(p) = u [xf{p)], 
S(x,p) = (l/3)pduf/dp — (1/3) J^° u x pdG J dpdx and 
G(x,p) — lim^oo g(t;,x,p) = S/w. Eq.( 20) describes the 



transition in the particle distribution between its asymptotic 
value g = G at £ — > oo and g — at £ — * — oo as a result of 
particle losses caused by the lack of resonant waves towards 
the subshock. The position of the transition front (£(x,p) = 
0) is determined by the condition of non-secular behaviour 
of the next order asymptotic expansion and must coincide 



with one of the characteristics of the operator on the l.h.s. of 
eq.@, i.e., 
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The reasonable choice of the concrete characteristic is based 
on the existence of the absolute maximum momentum p max 
beyond which there are neither particles nor waves. That 
means U{(p) = u [xf (p)] — uip/p max . Likewise, the func- 
tion x = Xf(p) is defined by Xf (p) = u^ 1 (uip/pmax)- 

3.0.1 External solution 

While having obtained the form and the position x = Xf (p) 
of the narrow front in the particle distribution g(x,p) we still 
need to calculate g to the right from the front where it decays 
with x. This would be the external solution G(x,p) intro- 
duced above. It is clear that 

max g(x,p) « G(xt,p) = G (p) 



(21) 



so that from eq.([18j) we have the following equation 
( \r i\ l r dG 
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The most important information about G(x,p) is contained 
in Go(p) for which from the last equation we obtain 
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where we have introduced v(p) by 
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Eq.(22) can be easily solved for Go 
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(where C is a constant). However, the function v depends on 
the solution itself. Nevertheless, it can be calculated prior to 
determining Go and therefore, this solution may be written 
in a closed form. In the case p ~ p max one obtains (the shape 
of the cut-off) 



Go ~ (Pmax - P) 



(25) 



In the rest of the a:,p-domain where X{(p) < x < L c and 
p is not close to p max , we may assume that the CR diffu- 
sion coefficient is close to its Bohm value since in contrast 
to the phase space region x rs X{(j>) at each given x,p there 
are waves generated along the entire characteristic of eq.(n3) 



passing through this point of the phase space and occupying 
an extended region of the CR precursor, Fig. 1 . We may adopt 
then the asymptotic high Mach number solution (Malkov and 
Drury 2001 [MD] and references therein) 

G(x,p) ~ Go (p) exp / udx I 

V k b J x , J 

Using the linear approximation for u{x) [MD] u — uq + 
mx/L' c with L' c = kb(P)/ui, where it is implied that the 
maximum contribution to the particle pressure comes from 
some momentum p ~ p, we can express v in the form of an 
error function integral 
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The further algebra simplifies in two limiting cases leading 
to the following asymptotic expressions for the particle dis- 
tribution 
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This result is valid forp > p^ = p max /R = p max 

whereas forp < p* one still has Go = C/p q ^~ 3 , where the 

detailed calculations of q(p), C and R may be found in MD. 

4 Conclusions 

Refraction to shorter wave lengths in a nonlinearly modi- 
fied flow results in a spectral break at p = p* = p ma x/R 
where R = Ui/uq 3> 1 is the nonlinear pre-compression 
of the flow. The spectrum rollover (p ma x/R < P < Pmax) 
is described by / oc p~ g e - v P'^with even faster decay at 
p ~ Pmax while an approximate power-law p~9(p) is still 
valid at p < p* . Since R itself is proportional to p* , the spec- 
tral break p* should grow slower than ^Jpmax- Due to the 
lack of waves at p ma x, it advances slower than in the Bohm 
case. 
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